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We propose a new exponential f{R) gravity model with f{R) — {R~ \c)e^^'^'^ and n > 3, A > 
l,c > to explain late-time acceleration of the universe. At the high curvature region, the model 

*vj ' behaves like the ACDM model. In the asymptotic future, it reaches a stable de-Sitter spacetime. It 

is a cosmologically viable model and can evade the local gravity constraints easily. This model share 

(— ^ , many features with other f{R) dark energy models like Hu-Sawicki model and Exponential gravity 

r\l ■ model. In it the dark energy equation of state is of an oscillating form and can cross phantom 

divide line ujde = — 1. In particular, in the parameter range 3<n<4, A~l, the model is most 
O „ distinguishable from other models. For instance, when n = 4, A = 1, the dark energy equation of 

state will cross —1 in the earlier future and has a stronger oscillating form than the other models, 
the dark energy density in asymptotical future is smaller than the one in the high curvature region. 

f^ \ This new model can evade the local gravity tests easily when n > 3 and A > 1. 

PACS numbers: 98.80.Hw, 04.80.Cc 
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^ ■ I. INTRODUCTION 

O . As we know, the standard big-bang cosmology based on radiation and matter dominated epochs can be well described 

^ ' within the framework of General Relativity [l|, |2|. The rapid development of observational cosmology starting from 

^ ' 1990s shows that the expansion of our universe in the present epoch is accelerating. Currently, the most successful 

I I model of cosmology that we have is the ACDM model. It is in well match with a wide variety of modern cosmological 

observations that have stunned the physicist community in the last decade. However, the ACDM model is not perfect 

C^ , in many aspects. First of all, the cosmological constant remains a mystery, can not be explained clearly in any 

known theory. If the cosmological constant originates from the vacuum energy in quantum field theory, as many 

people believe, its energy scale is too large to be compatible with the observed dark energy density Q. Moreover, 

[^ ' the observation indicates that the dark energy equation of the state may cross the phantom divide line w = — 1. 

^ ^ "^ 



This suggests that the cosmological constant [J-lyl may not be the only candidate for dark energy. There have been 
proposed many dynamical dark energy models to explain cosmic acceleration, ranging from quintessence, phantom, 
quintom to chaplygin gas models. For the nice reviews on dark energy, see [Tj-j^]. In dynamical dark energy models, 
p^ . one has to introduce at least one dynamical scalar to drive late-time acceleration, similar to the scalar driving the 
early-time inflation. 

An alternative scenario for dark energy is infrared(IR) modified gravity. Among many IR modified gravity model, 
f{R) gravity is of particular interest. One important feature in f{R) gravity is the intrinsic existence of an extra 
dynamical scalar degree of freedom, besides the massless graviton. Therefore it is possible to study both the early- 
?-j ■ time inflation and late-time acceleration of the universe in the framework of f{R) gravity, without introducing ad hoc 
scalar fields by hand. More interestingly, the effective equation of state could be smaller than —1 in f{R) dark energy 
models, indicating the scalar behaves like a phantom in the Jordan frame. On the other hand, the existence of scalar 
mode is not always pleasant. The fact that the dynamical scalar field may induce a long-range fifth force suggests 
that a viable f{R) gravity should satisfy the stringent constraints of local solar system test. 

Since the discovery of late-time acceleration of the universe in 1998, the f{R) theory have been extensively studied 
as the simplest modified gravity scenario to drive late-time acceleration. The model with a Lagrangian density 
f{R) = R — alR^ (a > 0,n > 0) was proposed for dark energv [l0l - [l3{ . However this model is plagued by matter 
instability |14l . ll5l and difficulty to satisfy local gravity constraints. Later on, researchers have proposed many viable 
models, seeing |16l423l ]. The Lagrangian of these models have a common form, adding a function of g{R/ Ri^) to the 
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Einstein-Hilbert term R. In the high curvature region, when R^ Rq, g{R/Ro) tends to be a constant and the model 
mimic the ACDM model, and should satisfy the local gravity constraints. However, in the late-time universe, the 
extra terms g{R/Ro) may play a significant role in the evolution. According to the dynamics of the theory, the dark 
energy equation of state could cross the phantom divide line, and tends to be —1 in the asymptotic future. For the 
nice reviews on f(R) theories, see |24l.l25|. 

In this paper, we propose a different f{R) dark energy model that do not contain a cosmological constant. In our 
model: 

/(i?) = (i?-Ac)e^(^)", (1) 

where A > 1, c > 0. The stability condition at the asymptotic future requires that n > 3. Among three parameters, c 
has the same dimension as Ricci scalar, A and n are dimensionless parameters. Different from the usual f{R) models, 
the Lagrangian of our model can not be separated into a R + g{R) form. Obviously, in the high curvature region, 
the exponential factor tends to be 1 and the model reduces to the ACDM model. Asymptotically, there is a stable 
de Sitter vacuum, with a different cosmological constant. We study the cosmological implications of this model. We 
discuss when and at what level it modifies the cosmological predictions, while evade the local tests of gravity. It turns 
out that even though when A ^ 1 or n ^ 3, the model becomes indistinguishable from the ACDM model, the model 
present distinguishable features from other models when A = 1 and 3 < n < 4. We compare our model with other 
two well-studied dark energy f{R) models, Hu-Sawicki modelflal and Exponential gravity modelfl^, |20,|23]- We find 
that all of them share some common qualitative features: crossing phantom divide line ujde = — 1, and dark energy 
equation of state being of an oscillating form, but they differ in details. 

This paper is organized as follows. In Section II, after briefly reviewing general f{R) theory, we introduce our model 
and discuss its cosmological implications. Via numerical analysis, we study the evolution of the FRW universe in our 
model and investigate the evolutions of the dark energy density and dark energy equation of state. In Section III, we 
analyze the local gravity constraints on our model. Finally, in Section VI, we present our conclusions and discussion. 

II. MODEL AND COSMOLOGICAL IMPLICATION 
A. Model 

The action of modified f{R) gravity with matter is 

S^ /d4xy^(ffl+L„), (2) 



2k^ 

where g is the determinant of the metric and R is the Ricci scalar curvature. Taking the variation of the action ([2]) 
with respect to gfj,,^, we have the equations of motion 

R^. - ^Rg^. = jT;Z + ;^{^5m.[/(^) - Rf] + (v^v. - 5^.0)/'}. (3) 

Here R^^i, is the Ricci tensor, /' — df /dR, V^ is the covariant derivative operator associated with the metric g^i,, and 
Dcf) = g^'^W f^W ^(f). And T™ is the matter stress-energy tensor which satisfies the continuity equation 

V^T™ = 0. (4) 

The trace of Eq. ^ gives 

30/' + f'R - 2/ = K^T, (5) 

where T is the trace of T^. 

The Einstein gravity, with a cosmological constant, corresponds to f{R) = R — 2 A, /' = 1 and R — 4A — k^T. 
Especially, in the matter dominated epoch, R ~ —k?T. In a general modified gravity, /' could be considered as a new 
scalar degree of freedom, (f) oc /'. The trace equation ([S]) determines the dynamics of this scalar field (}>: 

□/' = ^, (6) 



with the effective potential 



9K 
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In the late-time universe, p <?C 1 can be neglected, the equation dVe/f/df 
d^Veff/df^ > gives 



and 



2f{R) - Rf'iR) = 

'Veff _K.f' 



Of- 



3 7" 



R) >0. 



(7) 
and the stability condition 

(8) 
(9) 



If in a f{R) gravity the solution of Eq.® gives a positive scalar curvature and satisfies the stability condition ^, 
then the universe will enter into a stable de-Sitter phase in the asymptotic future. 
To be cosmologically viable, a f{R) gravity model should satisfy a few requirements: 

• In the high red-shift regime, f{R) — > i? — 2A, and mimic the ACDM model which is well tested by the CMB, 
supernova and other experiments; 

• In the asymptotic future, the model will have a stable vacuum, which is usually a de Sitter spacetime; 

• It should satisfy the local gravity constraints in the high curvature region. 

In the cosmologically viable f(R) models in the literature, they usually take a form as R + g{R/Ro), where g{R/Ro) 
tends to be a constant when R/Rq — >■ oo. We select two well-studied ones[la,l23 to compare with our model: 
Model f{R) Parameters 

(i) Hu-Sawicki R - ^^^^SIf^+i" ^i' ^2, P(> 0),i?Hs(> 0) 

(ii) Exponential R - /SRe (l - e-^/^^) ^, Re 

(iii) Our model (R - Xc)e^^Ti'>" A>1, ri>3, oO 



C2{R/Rhs)'' + 1 

R^l3RE{l-e-^/^^) 

In the matter-dominated epoch, the curvature is large and the exponential factor in our action could be safely set 
to 1 and the model reduces to ACDM model with A = ^. 

To study the asymptotic behavior of our model, we need to solve Eq.®, taking into account of the stability 
condition (j9]) . From the explicit form ([T|) , we calculate the first and the second derivative of / 
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If A is large enough or n » 3, the de-Sitter solution of Eq.® is 

Rds ^ 2Ac, 



(10) 
(11) 

(12) 



which is the same as the solution in the ACDM. For general values of A and n, Rds J^ 2Ac. We will discuss this issue 
via numerical analysis in the next section. 
On the other hand, the stability condition 



9/'2 3 "'^f'Rds 
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-Rds 






3 -yXnc-[n{l-^)-{l + ^)] 
is easily satisfied when n > (1 + -^)/(l — 7^) — 3. If we define 

f"{R)R 



(13) 



(14) 
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FIG. 1: The parameter m{Rds) as a function of A with n = 3, 4, 5, 6 at the de-Sitter universe. 



we find that the stability condition becomes simply < ?7i < 1 , noticing that the Compton wavelength of the extra 
scalar mode is defined as 






(15) 



/37"~ ^/3m/R. 



In order to analyze the stability precisely, we have plotted in Fig. 1 the value of parameter m in the asymptotic 
de-Sitter phase for different values of n. When n > 3, < m{Rds) < 1 such that the stability condition is well 
satisfied. However, when n < 3, the model is unstable. Another feature learned from the Fig.l is that when A > 4, 
the value of m is very tiny. 

We will work in the FRW space-time. In this case, Eq. ([3|) gives the modified Friedman equations: 



1 



3f'H^ = K^p.^ + -{f'R-f)-3Hf', 



(16) 

(17) 



in which we have neglected the radiation component because we only focus on the cosmology from matter-dominated 
epoch to the asymptotic future in this paper. The dot denotes time derivative d/dt. The matter density satisfies the 
conservation law 



and can be written as p„ 



Poa 



Prn + iH{p^) = 0, (18) 

■^. The dark energy density and pressure can be defined respectively as: 



1 



^'Pde = f' + 2Hf'-^{fR-f)-{2H + 3H^){l-n. 



(19) 
(20) 



B. Dynamical evolution 

According to the work of Linder|20ll. with some modifications, we define 

H"^ _, Ac 



Xh = ^-a ^-TTT' (21) 

mg DTTlg 

R . \c 

XR - -^-2,a-^ -2^.-l2xH, 22 

mi mf. 



'0 

where the parameter TOq is defined as k^pq/3. The matter density parameter is defined as 



"" = ^1 — (23) 

Using Eq. ((T6|) . we get 

3/'i/2 = «:V,„ + i(/'i? - /) - 3i7V"^, (24) 

where A^ = Ina. Using R = 12H^ + 3d{H^)/dN, the modified Friedman equations become two first-order equations 

dxH XR 



dN 3 



(25) 



Here, R and _ff can be expressed in terms of the parameters xh and xr. The dark energy equation of state is 
defined as 



-2H - 3H^ _ 1 XR 

3H^ — K^Pm 9 XH + Ac/6mg 



LJde= ^^^., =-l-7T ^__ , x:/^^2 - (27) 



Now let us study the expansion history of the universe in our model. At the high curvature region dominated by 
matter, according to the trace equation ([5]), the new scalar /' is always near the minimum of the effective potential, 
the Ricci scalar is approximatively. So, we give the expression of Ricci scalar 

i?~K2pm + 2Ac + e'(i?(c/i?)"). (28) 

On the other hand, f{R) -^ R — Xc mimic R — 2A. Using Eq. ([T6|) . and taking the initial value of /' as zero, the 
Hubble parameter can be written as 

H^^^ + ^+0{R{c/Rr). (29) 

3 6 

In the early universe, Xh and xr^ are the same order as R{c/R)^ and can be neglected. Consequently, we may take 
the initial condition xh = xr, = at z=13. 

In the asymptotic future, the universe will go into a de-Sitter epoch as expected. The stationary point of the 
equation of motion is given as 

XR = 0, (30) 

\{f'R~f)-if'H'^^Q, (31) 

which is equivalent to Eq. (jS]) . The fact that xr evolves from zero to zero suggests that the equation of state of dark 
energy starts from —I at the matter-dominated epoch and will end with the same value —1 in the future. With the 
matter density becoming smaller and smaller, the Ricci scalar would tend to be a constant Rds- 
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FIG. 2; Density parameters Q,„ and Oa in our model with the parameters A = 1, n=4. 




FIG. 3: Dark energy density p^e, normalized by K^po/3. 



C. Numerical analysis 



From Eq. ([23)) . we plot the density parameters flm and O^e as a function of Ina in Fig. 2. Note that, the current 
matter density parameter fl^ — 0-3 is in accordance with the observation. 

In Fig. 3 we show the evolution of the dark energy density in our model. In the high red-shift region A^ < — 1, 
the dark energy density stays around the value Ac/2. In the region — 1 < A^ < 1, it first rise up to a maximum, and 
then falls down and oscillates around another constant value asymptotically. Note that the dark energy density tends 
asymptotically to different values with different parameters A. A smaller A leads to a smaller asymptotic value of dark 
energy density in the future. When A is large enough, the k? pde would tend to be Ac/2, as the same as the value in 
the matter-dominated epoch. 

In Fig. 4, we show the evolution of the Ricci scalar curvature for different values of A. We see that R evolves from 
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FIG. 4: The Ricci scalar curvature, normalized by k po/3 evolves from high curvature to a constant. 
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FIG. 5: The dark energy equation of state of for different models. 



a high curvature and quickly reaches its asymptotic value Rds- Note that a smaller A will lead to a smaller value of 

Rds- 

In Fig. 5, we show the evolution of the equation of state of dark energy for different models and parameters. From 
the dark energy conservation equation, we know Ude — —1 — Pde/{iHpde)- We can see the relationship between w^e 
and pde from Fig. 3 and Fig. 5. At the beginning, ujde stays near —1 in high red-shift region N < —1. It first falls 
down to a minimum and then climbs up to a maximum. It falls down again and begin oscillating around —1 with 
smaller amplitude. It finally settle down to —1. We compare the models of Linder and Hu-Saweicki with ours. We 
find that in all models, the evolution of dark energy equation of state is quite similar qualitatively. However, the 
details of the evolution are different. In particular, when n < 4, the oscillation behavior of our model looks different. 
For example, in the case n = 4, in our model ujde cross phantom divide line —1 earlier and oscillates in the region 
— I < In a < 1. Fig. 5 indicates also that a larger value of the parameter n will reduce the deviation from —1, 
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FIG. 6: The parameter m, as a function of Ina. 
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FIG. 7: The parameter m as a function of Ina with parameter n ~ 3. 



seeing the red lines with parameters n = 4,5,6,A = l. So does a larger value of A, seeing the lines with parameters 
A = 2,1.5,l,ri = 4. The fact that larger values of n and A will suppress the deviation of Wde from —1 is easy to 
understand: the larger n or A, the closer our model to the ACDM model. 

In Fig. 6 and Fig. 7, we show the evolution of the parameter m which decides the Compton wavelength Ac = 
■^3m/R. The parameter m in the models of Hu-Saweicki and Linder starts from a small quantity, then climbs up to 
a constant. In our model, when n < 4, the oscillation of the parameter m in the evolution is obvious. This difference 
accounts for the fact that the dark energy equation of state in our model has a different oscillating form from the ones 
in the other two models. When n = 5, 6 the shape of the evolution lines is similar to the ones in the other two models. 
In general, a larger A will lead to a smaller m. When n — ;► 3, the terms of (c/R)" and (c/i?)""*""'^ in the expression of /" 
cancel out each other, the terms of (c/R)'^'^ and (c/i?)^"+^ become important. Therefore the parameter m becomes 
even smaller asymptotically, seeing Fig. 7. 
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FIG. 8: The square of Hubble parameter in our model comparing with ACDM, and its evolution. 

In Fig. 8, we show the evolution of the Hubble parameter in our model, comparing with the ACDM model. We see, 
in the high red-shift region, the ratio tends to be 1, which means the f{R) is close to the ACDM. However in the low 
red-shift region, it will be a little larger than that in ACDM. And in the future — 1 < z < 0, the ratio will be less than 
1. Note that, a smaller value of the parameter A suggests a larger deviation from ACDM in the late-time universe. 
As a result, the current age of the universe in our model is smaller than the one in the ACDM model. Moreover, the 
cosmological distance is also be affected. Thus we can use the observational data to constrain the parameters. For 
n = 4, the deviation from the ACDM model appears at z = 3. 

III. THE LOCAL GRAVITY TEST 



In this section, we discuss the compatibility of our model with the local gravity test. As a f{R) gravity theory is 
equivalent to a scalar-tensor theory, the extra scalar mode may mediate a long range attractive fifth force and thus 
violate solar system constraints. However, it has been suggested in (la. [26l - l3l| . the /(i?) models can be consistent 
with the local gravity constraints with the help of the chameleon mechanism. In the chameleon mechanism, the scalar 
behaves differently in different environments. This is achievable in f{R) gravity as the scalar potential get modified 
by the scalar coupling to the matter density. 

It is more convenient to work in the Einstein frame, which is related to the original Jordan frame by a Weyl scaling 
5^1/ = f'Stiv- The extra scalar field (p is defined as 2/3k0 = In/' with f3 = l/v6- The action in the Einstein frame is 



S 



EF 




1 



^^..^^(p^,(p-Vm+ / d^XLm{l}f..f-\^M) 



where 



F(0)- 



I'R-I 



2^2/ 



2 f /2 



(32) 



(33) 



If we consider the background geometry as a Minkowski spacetime, then the equation of motion for the scalar </) in 
the Einstein frame is just 



dK 



e// 







d/P' r dr d(f> 
where r is the distance from the center of the object, and Veff is the effective potential, 

Veff=V{cf>)+p,ne-''^t 



(34) 



(35) 
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It is this form of the potential that makes the chameleon mechanism feasible. 

Let us consider a spherically symmetric object with a radius r^, mass M,,, a matter density p — pa when r Kr^i 
and p = pb when r > r^. The effective potential has a different shape with a different surrounding matter density. It 
has minima at (pA and (J)b'- 

V,^ {(J^a) - PAPne-^^^^ = 0, (36) 

V,^ (0b) - Ps/3Ke-^«^^ = 0, (37) 

where the scalar has different mass m^ — d'^Veff{(j)A)/d(p'^ and m^ = d'^Veff{4>B)ld4P' respectively. For our model, 
the first derivative of the effective potential is 

dVeff _dVdR ^_ ^_^,^ 



^9i?9^~^"^"^" 



13k 



(2/ - I'R) 



2 f /3 



K^f 



Pme 



-,3k0 



(38) 



; 2f-f'R 



■m , 



where we have used /' = e'^^'^'^ ~ 1 and /" = 2/3 /te^^"*^^ because that the minimum of the potential is very near 
in the high curvature region Rq ~ k? Pm S> 1. According to the definition we have, 



^ '" ■ \2^(' ^ \n+l 



2/3k0™„ = In/' ~ -A(n - 1)( — )" + A2n( — )"+\ (39) 

The second derivative of the potential gives the scalar mass 

If the object satisfies the thin-shell condition k((/)s — (?!)a)/6/3<I'c ^ 1; where $c is the Newton potential, then the 
scalar has an exterior solution 

If the object is the Sun, we have pA — Ig/cm^ and pb — \Q~'^^g/cm^ which is the galaxy matter density. The 
cosmological density is po = \Q~^^g/cw? . Thus Ra — k?Pa is much larger than Rb, and according to Eq. ((38|) . we 
have 

2/3K0B^-A-"+i(n-l)(^)" 

k-^Pb 

~ -A-"+i(n - l)(4^)"10-5" (42) 

:^_A-"+i(n-l)(^)"10~5", 

where we have used the fact that Ac — > 2 A. Obviously (t)B is much larger than 4)a which can be neglected in Eq. (|4ip . 
From Eq. (HH)) . we have 



3m(_RB) 

"a V ^n 



10"^^^^ /2fU J ^43) 



10 



-16 



where we have used $c = GMc/tc = k^pa^I/Q — 10 ^. We see that the mass of the scalar is very light out of the 
object, so the exponential e"^^^''"'''^^ in Eq. (PT|) could be set to 1. 
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Let us transform back to the Jordan frame. Under the inverse transformation, g^i, = e~'^^'^'^'g^j,^, r = e^^'^'^r, the 
metric in the Jordan frame is 

ds^ = e-2/3«0dS^ - -[1 - 2A{r)]df + [1 + 2B]dr^ + r^d^"^. (44) 

under the condition j3K(f> <C 1. Then we have the following relations 

Air) ~ A(r) +/?«</, ~ —^[1 + -^i- - 1)], (45) 

Bir) ^ Bi?) + Pk?^ ^ ^^(1 + pA. (46) 

dr r 9r 



The tightest experimental bound on the PNP parameters is given by|7— l|<2.3xl0 ^ [24l l35l - l37| . If we take the 
distance r = rc^ then we can constrain the parameter in our model 

|/3k0b| = iA-"+i(n - l)(^)"l0-5" < 2.3 x 10"". (47) 

As n > 3 and A > 1 in our model, this bound can be satisfied easily. 

On the other hand, the experiment from the violation of equivalence principle gives a slightly more stringent 

bound m US miii, 

1^1 < 8.8x10-7/3, (48) 

which can be translated into 

|/3k0b| <2.1 X 10"^^ (49) 

It could be evaded in our model without trouble. 

IV. CONCLUSIONS AND DISCUSSION 

In this paper, we have proposed a new viable fJR) dark energy model. It is of an exponential form, but is 
different from the Exponential gravity proposed in [20| . We focus on the cosmological evolution starting from matter- 
dominated epoch to the asymptotic future. In the matter-dominated epoch, our model reduces to the ACDM model 
with a positive cosmological constant. In the asymptotic future, the universe will settle down to a stable de-Sitter 
phase. However, between two epoches, the evolution of the universe in our model could be very different from the 
ACDM model, if we choose the parameters appropriately. 

The solar system constraints of f{R) gravity place weak bounds on our model. Due to the exponential suppression, 
our model has little deviation from GR in the high curvature region. It can evade local gravity constraints easily. 

It turns out that in the parameter range 3<n<4, A^l, the model is most distinguishable from other models. In 
our model, the dark energy equation of state cross the phantom divide line many times, before it asymptotically settle 
down to a constant —1. Comparing with the other f{R) models, the dark energy equation of state in our model has 
strong oscillations, and will cross the phantom divide line in the earlier future. This prediction can be tested in the 
future observation. In our model, the Hubble parameter shows the deviation from the ACDM model at the redshift 
z « 3, suggesting that the current age of the universe and the theoretical distance of the cosmological objects may be 
smaller in our model. It would be very interesting to further constrain the parameters in our model from cosmological 
observations. 
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